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We present results and predictions for the nuclear modification of the differential cross sections 
for inclusive light hadron and prompt photon production in minimum bias d+Au collisions at y/s = 
200 GeV and minimum bias p+Pb collisions at y/s = 5 TeV at RHIC and LHC, respectively. Our 
calculations combine the leading order perturbative QCD formalism with cold nuclear matter effects 
that arise from the elastic, inelastic and coherent multiple scattering of partons in large nuclei. We 
find that a theoretical approach that includes the isospin effect, Cronin effect, cold nuclear matter 
CI energy loss and dynamical shadowing can describe the RHIC d+Au data rather well. The LHC 

p+Pb predictions will soon be confronted by new experimental results to help clarify the magnitude 
D ' and origin of cold nuclear matter effects and facilitate precision dense QCD matter tomography. 
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I. INTRODUCTION 
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Medium-induced modification of high transverse momentum (pt) particle production in nucleus collisions (p+A, 
A+A) relative to the naive binary collision-scaled proton-proton (p+p) baseline expectation is a powerful probe of 

l— '■ the properties of dense QCD matter. In particular, the nuclear modification in d+Au collisions at the RHIC and the 
forthcoming p+Pb results at the LHC can provide valuable information on the elastic, inelastic and coherent multiple 
parton scattering processes inside a large nucleus and are vital testing grounds for novel nontrivial QCD dynamics [l[ . 
The manifestation of such nontrivial QCD dynamics in p+A reactions is usually referred to as cold nuclear matter 

£T) . (CNM) effects, as opposed to the effects of parton and particle formation and propagation in the ambiance of the 

CZ^ quark-gluon plasma created in high energy A+A collisions [2] . 

^^O , There have been several different approaches to studying CNM effects. One of them is based on the leading- twist 
^\ ■ perturbative QCD factorization and attributes all these effects to universal nuclear parton distribution functions 
(nPDFs), which are the only ingredients different from the case of p+p collisions [3|-l5|- Through a global fitting 
procedure to fix the parametrization of nPDFs, this approach can describe part of the RHIC d+Au data reasonably 
well and has made predictions for p+Pb collisions Another approach is the so-called Color Glass Condensate 

J> . (CGC) approach [fj-[lj]. It focuses on non- linear corrections to QCD evolution equations in very dense gluonic systems 
and is only applicable in the very small- a; region - the so-called classical gluon saturation regime. Various predictions 
for the p+Pb run have also been made based on CGC approach, see Refs. p^4Tl |. There are also calculations and 
■ predictions for the nuclear modification factor in Monte Carlo models such as HIJING [l5| ■ 

In this manuscript we follow an approach different from the above mentioned. It is based on perturbative QCD 
factorization and CNM effects are implemented separately within the formalism. The advantage of this approach is 
that all CNM effects have clear physical origin, mostly centered around the idea of multiple parton scattering 16!]. 
Thus, they can be calculated and their implementations is well documented in the literature. They are also process- 
dependent. In this manuscript we discuss the isospin effect, Cronin effect, cold nuclear matter energy loss, and 
dynamical shadowing. The isospin effect is easily incorporated by replacing the nPDFs by the Z, A — Z-weighted 
average of proton and neutron PDFs [l7], EH ■ Theoretical interpretations of the Cronin effect are reviewed in [l9| 
and our implementation is via initial-state parton transverse momentum broadening. Initial-state cold nuclear matter 
parton energy loss was calculated in pol l2l| . Nuclear shadowing can be implemented as coherent power corrections 
to the differential particle production cross section [22], [H| ■ A perturbative QCD formalism with these effects incor- 
porated can be used to understand the qualitative and in most cases quantitative features of the nuclear modification 
in p+A (or d+A) collisions observed at CERN SPS and RHIC experiments [231 - 4261 ] . In anticipation of the new LHC 
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p+Pb results we present predictions in the framework of this approach to incorporating CNM effects. Specifically, 
we focus on inclusive particle production in d+Au and p+Pb collisions and show results for the nuclear modification 
factor of neutral pions, charged hadrons and prompt photons. 

This paper is organized as follows. In Sec. II, we present our approach to studying CNM effects. We first review 
the basic pcrturbativc QCD formalism for single inclusive particle production in p+p collisions. We then discuss how 
various CNM effects are incorporated in this formalism. In Sec. Ill, we implement these CNM effects and show that 
the formalism can describe particle production in d+Au collisions at RHIC energies. We then present our theoretical 
model predictions for the upcoming p+Pb collision results at the LHC. We summarize our paper in Sec. IV. 



II. SINGLE INCLUSIVE PARTICLE PRODUCTION 



In this section we review the basic perturbative QCD formalism for single inclusive particle production in p+p 
collisions. In particular, we focus on light hadron and prompt photon production. We then discuss the incorporation 
of various cold nuclear matter effects within this formalism. Specifically, we comment on the isospin effect, Cronin 
effect, cold nuclear matter energy loss and dynamical shadowing. 



A. Basic pQCD formalism 



To leading order in the framework of factorized perturbative QCD, single inclusive hadron production in p+p 
collisions, p(pi) + p(pi) — > h(ph) + X, can be written as follows (27| : 



da 



= K—^2 [— d 2 k aT f a /N{x a ,kl T ) [— d 2 k bT f b/N {x bl kl T ) 

5 , J X a J X b 



dyd 2 p T 

dz c 

D h / c (z c )Hab^c{S 7 t, u)S(s + t + u), (1) 

z c 

where y and px are the rapidity and transverse momentum of the produced hadron and ^2 a b c runs over all parton 
flavors. In Eq. (p} s = (pi +P2) 2 , D h / C {z c ) is the fragmentation function (FF) of parton c into hadron h, H ab ^ c (s, i, u) 
are the hard-scattering coefficient functions with s,i, u being the usual partonic Mandelstam variables [13, [Hj]. We 
also include a phenomenological if-factor to account for higher order QCD contributions f a ,b/N(x, &f0 are the 

parton distribution functions with longitudinal momentum fraction x and transverse momentum component fcy- We 
have included this fey-dependence in order to incorporate the Cronin effect in p+A collisions. We assume a Gaussian 
form in this variable [I3,[27|], 

fa/N(x a ,k 2 aT ) = f a/N ( Xa )-±^ e - k ^/^\ (2) 

where f a /N( x a) are the usual collinear PDFs in a nucleon. 

On the other hand, prom pt p hoton production in p+p collisions has two components - the so-called direct and 
fragmentation contributions |27l . |30| : 

da _ da diT - + da frag - 
dyd 2 pT dyd 2 px dyd 2 pr 

The fragmentation contribution is given by: 

dcr fra K . „,2 



dyd 2 pr 



= K ~ "52 f ' —d 2 k a T fa/N(Xa,kl T ) [ — d 2 k bT f b/N (x b , k\ T ) 
S u J X a J X b 



a,b,c ' 

dz c 

D 1 / c {z c )H a b^c{s 1 t, u)8{s + t + u). (4) 



In other words, it is exactly the same as Eq. (JT|) if one replaces the parton-to-hadron FF Z?^/ C (z c ) by the parton-to- 
photon FF D 1 / C {z c ). On the other hand, the direct contribution can be written as: 

da — = K a oma s ^ I d 2 k aT f a / N {x ai kl T ) f —d 2 k bT f b/N (x b , k bT ) H ab ^ 7 (s,i,u)S (s + i + u) , (5) 
dyd pt s j Xa J x^ 

a.b 
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where _ff a h_j. 7 are the well-known partonic hard-scattering functions for direct pho ton production [27J, |3l|, [32| . 

In the numerical calculations we choose (k^} pp = 1.8 GeV 2 in p+p collisions |24l [25j| . We use CTEQ6L1 PDFs [33j . 
fDSS parametrization for parton-to-hadron FFs [HI and GRV parametrization for parton-to-photon FFs (35[ . We fix 
the factorization and renormalization scales to the transverse momentum of the produced particle, fif — fi r = pt- 
We find that a 0(1) if- factor can give a good description of hadron production in both RHIC and LHC energies, 
see Fig. [TJ For prompt photons the if- factor values are slightly larger and also -^/i-dependent but, importantly, 
the shape is well-described. In the left panel of Fig. [TJ we compare our calculation to RHIC photon and ir° data 
at y/s = 200 GeV [H, [37[ ■ We also consider LHC photon data at y/s — 7 TeV [38| and charged hadron data at 
\/s = 2.76 TeV [39j . shown in the right panel. The leading order pQCD formalism provides a reasonable baseline 
description of particle production in p+p collisions. The precise value of the if -factor is not important for the study of 
the nuclear modification of inclusive particle cross sections in p+A collisions (it cancels in ratios) . We will incorporate 
nuclear effects relevant to p+A reactions in this formalism in the next subsection. 



O 



T3 

ffl 



io! 

10" 



10' 



Qh 
T3 



10 
10 



10 



10 ' 
10 4 
10 

-2 

10 



200 GeV 

PHENIX photons 
Theory K=l.f 




Vs=200 GeV 
PHENIX 7i°s 
Theory K=l 



10 



15 



20 



25 30 
P T (GeV) 



a 

H 

-3 



10 



10 



> 
a 



-d 

W 



10 



-5 - 



10 



10 




7 % 



10 



-12 - 



Vs=7 TeV 

o CMS photons 
- Theory K= 1.5 



Vs=2.76 TeV 
CMS charged hadrons 
- - Theory K= 1 



V 



-■a 



50 100 150 200 250 300 350 400 

P T (GeV) 



FIG. 1. Left panel: our calculation at \fs = 200 GeV and rapidity y — is compared to RHIC photon data (top) [3(| and 7r° 
data (bottom) [33]. We choose K — 1 for 7r° and K — 1.8 for photon production. Right panel: comparison to LHC photon 
data at y/s = 7 TeV and \y\ < 0.9 38] and charged hadron data at ^fs = 2.76 TeV and \y\ < 1 f39j. We choose K = 1.5 for 
photon and if = 1 for charged hadron production. 



B. Cold nuclear matter effects 



The p+A (e.g. d+Au or p+Pb) nuclear modification factor R p a is usually defined as: 

p = 1 da vA I da pp 
pA (JVcou) dydPpr/ dydPpr' 1 ' 

where (N co u) is the average number of binary collisions. The deviation of R p a from unity reveals the presence of CNM 
effects in p+A collisions. A variety of CNM effects can affect particle production. In this manuscript we discuss the 
ones that have been theoretically evaluated as arising from the elastic, inelastic and coherent scattering of partons in 
large nuclei We also account for the proton and neutron composition of the interacting nucleus. In particular, we 
implement the isospin effect, Cronin effect, cold nuclear matter energy loss and dynamical shadowing. These effects 
have been well documented in the literature. We give a brief overview of their implementation in preparation for LHC 
predictions. 

Isospin effect. Isos pin effect can be easily accounted for on average in the nPDFs for a nucleus with atomic mass 
A and Z protons via [17], 0: 



fa/A{x) = ^fa/p{x) + U - fa/n{x). 



(7) 



In Eq. f a / p (x) and f a / n ( x ) are the PDFs inside a proton and neutron, respectively. The PDFs in the neutron are 
related to the PDFs in the proton via isospin symmetry. 
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Cronin effect. Theoretical approaches to understanding the Cronin effect are well documented in jl9|. It can be 
modeled via multiple initial-state scatterings of the partons in cold nuclei and the corresponding induced parton 
transverse momentum broadening [Z(l EJ- In particular, if the parton distribution function fb/A( x b, kl T ) has a 
normalized Gaussian form, the random elastic scattering induces further fey-broadening in the nucleus (l7l l29j : 

Kt)vA = {%,t)pp + (^)s- (8) 

Here kb y r is the transverse momentum component for the parton prior to the hard scattering, ( = ln(l + 5p^) (2^.[25|. 
and we choose 5 = 0.14 GeV~ 2 , /i 2 = 0.12 GeV 2 , \ g = Cf/Ca^q = 1 fm. These parameters can describe reasonably 
well RHIC data. 

Cold nuclear matter energy loss. As the parton from the proton undergoes multiple scattering in the nucleus 
before the hard collisions, it can lose energy due to medium-induced gluon bremsstrahlung. This effect can be easily 
implemented as a momentum fraction shift in the PDFs 



fg/pM -> f q /p ^ 1 _ £cff J ' f B /p( x *) ~> fg/P ^ i_ eoff ) ■ (9) 

Ideally, Eq. (JSJ) should include a convolution over the probability distribution of cold nuclear matter energy loss (43| . 
However, concurrent implementation of such distribution together with the Cronin effect and coherent power correc- 
tions is computationally very demanding. The main effect of the fluctuations due to multiple gluon emission is an 
effective reduced fractional energy loss e c s. relative to the mean value (e) = ^§ L )- Here, the sum runs over all 
medium-induced gluons. In this work we follow Ref. [IH and use e c g. = 0.7(e). The average cold nuclear matter 
energy loss is obtained by integrating the initial-state medium-induced bremsstrahlung spectrum, first derived in [2(| ■ 
It also depends on the typical transverse momentum transfer squared /i 2 per interaction between the parton and the 
medium and the gluon mean free path X g . For this reason, the parameters are constrained to be the same as in the 
implementation of the Cronin effect, /i 2 — 0.12 GeV 2 and X g = 1 fm in our comparison to RHIC data in the next 
section. This calculation of initial-state cold nuclear matter energy loss has been shown to give a good description of 



the nuclear modification of Drell-Yan production in fixed target experiments 43 1. 

Dynamical shadowing. Power-suppressed resummed coherent final-sate scattering of the struck partons leads to 
shadowing effects (suppression of the cross section in the small- x region) (23[. The effect can be interpreted as a 
generation of dynamical parton mass in the background gluon field of the nucleus [44| . It is included via: 

x b ^x b (l + c / 2 ^.' 1 ^ , (10) 

where Xb is the parton momentum fraction inside the target nucleus, Ca — Cf{Ca) if the parton d = q(g) in the 
partonic scattering ab — > cd. £ 2 represents a characteristic scale of the multiple scattering per nucleon. In the RHIC 
energy range y/s = 200 GeV, £ 2 = Cf/Ca^ = 0.12 GeV 2 [22|,[23[ can give a good description for nuclear modification 
in d+ Au collisions for both single hadron and di-hadron production [23], [26| . 



III. THEORETICAL MODEL PREDICTIONS IN P+A COLLISIONS 

In this section we compare our theoretical model predictions to existing light hadron and direct photon data in 
minimum bias d+Au reactions at RHIC at y/s = 200 GeV. We also present predictions for these inclusive final states 
in minimum bias p+Pb collisions at the LHC at y/s = 5 TeV. It should be understood that in nucleus collisions we 
quote the center-of-mass energy per nucleon pair. 



A. Nuclear modification factor in d+Au reactions at RHIC 



We have incorporated the cold nuclear matter effects discussed in the previous section into the leading order pQCD 
photon and hadron production formalism, Eqs. ([1} and , . Such theoretical approach is able to give a good 
description of single inclusive particle production in d+Au collisions at RHIC energies. While such calculations have 
been presented in the past, in this manuscript we extend our results to direct photons at forward rapidity. In Fig. [2] 
we plot the minimum bias nuclear modification factor RdAu as a function of transverse momentum for both photon 
(left panel) and ir° (right panel) production at y/s = 200 GeV. The photon data is from the PHENIX collaboration 
at rapidity y = [36|. For n°s we have included both mid-rapidity y = [45| and forward rapidity 3 < y < 3.8 [46[ 
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FIG. 2. The nuclear modification factor RdAu is plotted as a function of transverse momentum pr for photons (left panel) 
and 7r°s (right panel) at RHIC energy of yfs = 200 GeV. Empty and filled circles are low-mss virtual photon and re al p hoton 
experimental data at mid-rapidity from PHENIX [3(|. Filled squares are the minimum bias mid-rapidity ir° data [45(|, filled 
triangles are the forward rapidity y = 4 STAR tt° data [43] , and empty squares is the forward rapidity 3 < y < 3.8 PHENIX 
7r°s [4y]. We have also presented predictions for rapidity y — 1.5 (dashed curves) and y — 3 (dotted curves) for both photon 
and hadron production. 



PHENIX data, and forward rapidity y — 4 STAR data 47]. We conclude that implementation of cold nuclear 
matter effects based on the physics of multiple parton scattering in dense QCD matter can describe both photon and 
hadron data quite well at mid and forward rapidities. With the future RHIC upgrade program in mind, we also show 
predictions for the nuclear modification factor RdAu for photon and 7r° production in different kinematic regions: the 
dashed curves are for y = 1.5 and the dotted curves are for y — 3. 



B. Predictions for the p+Pb run at the LHC 



We are now in good position to present predictions for minimum bias p+Pb reactions at the LHC energy of 
y/s = 5 TeV. These results will soon be confronted by new experimental data. The main parameters in our approach 
include: the Cronin momentum broadening parameters /i 2 , X q<g , which also determine the cold nuclear matter energy 
loss, and the dynamical shadowing parameter £ 2 . As emphasized in [31(, ^ 2 /A 9i9 and £ 2 represent the strength of 
the multiple scattering between the incoming and outgoing partons and the target nucleus, thus proportional to the 



number of the soft gluons in the nuclear medium. Following 31 [, beyond the possibility of fixed £ 2 , which we take as 



the lower limit for this parameter, we study a scenario where £ oc ai n , the inelastic nucleon-nucleon scattering cross 
section. Taking into account that Oj„ = 42 mb [48| at RHIC energy of y/s — 200 GeV and crj„ = 70 mb [H| at the 
LHC energy y/s = 5 TeV, we consider a potential increase in the £ 2 parameter by 67%. Specifically, for the upper 
limit of this parameter we take £ 2 = 0.20 GeV 2 in our prediction for the LHC p+Pb run at y/s — 5 TeV. We will also 
consider a correlated enhancement in /i 2 /A s . 

In Fig. [3] we present our theoretical model predictions for the nuclear modification factor R p pb as a function of px 
for neutral pion (left panel) and charged hadron (right panel) production in minimum bias p+Pb collisions at the 
LHC energy of yfs = 5 TeV. Rapidities y = (top), y — 2 (middle), and y = 4 (bottom) were considered. The upper 
edge of the bands corresponds to the RHIC parameters (a* 2 ,£ 2 ) = (0.12,0.12) GeV 2 . The lower edge corresponds to 
a potential enhancement of these parameters as discussed above, (// 2 ,£ 2 ) = (0.2,0.2) GeV 2 , to be tested against the 
experimental data. In Fig. 01 we present similar results for prompt photon production. 

The general features of these nuclear modification factors are very similar to the ones observed at RHIC energies. 
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FIG. 3. Predictions for the nuclear modification factor R p pb as a function of transverse momentum pr for tv° (left panel) and 
charged hadron (right panel) production in minimum bias p+Pb collisions at the LHC energy of yfs — 5 TeV. Results for three 
rapidities y = 0, y = 2, and y = 4 are shown. 



At mid-rapidity y — 0, there is a very small "Cronin peak" in the low px region. The peak is very close to unity 
and not as pronounced as the one observed in low energy fixed target experiments. This is because the dynamical 
shadowing becomes important and strongly suppress the particle production in this region. Furthermore, initial-state 
energy loss is larger due to the bigger contribution of diagrams with incoming gluons to the cross section. At high 
Pt we still have ~ 15% suppression, which is also due to the cold nuclear matter energy loss. In general, at forward 
rapidity all CNM effects are amplified due to the larger values of the momentum fraction x a from the incoming proton 
(relevant to cold nuclear matter energy loss) , the smaller values of the momentum fraction Xb of the incoming nucleus 
(relevant to dynamical shadowing) and the steeper falling spectra (relevant to the Cronin effect). As a result, at 
forward rapidity and at low px the dynamic shadowing effect can be most important and lead to stronger suppression 
of inclusive particle production while at high px the suppression is a combined effect of cold nuclear matter energy loss 
and the Cronin effect. Finally, the behavior of the R p pb for 7r°s, charged hadrons and direct photons is qualitatively 
the same and the quantitative differences are minor. 

It is interesting to point out that our predictions for inclusive particle suppression at forward rapidities are for 
observable effects that are stronger than the effects found in calculations using EPS09 nuclear PDFs. On the other 
hand, they lead to weaker suppression than the CGC predictions from a hybrid formalism in rc-BK-MC approach 0, 
14]. In other words, our predictions fall between EPS09 and CGC based approaches. 



IV. SUMMARY 



We studied the nuclear modification of high transverse momentum particle production in minimum bias d+Au and 
p+Pb collisions at RHIC and LHC, respectively. Such nuclear modification is manifestation of nontrivial and not 
yet fully understood QCD dynamics in large nuclei. It is, therefore, critical to further our understanding of high 
energy nuclear reactions through concurrent theoretical advances and comparison of model predictions to precise 
experimental data. With this motivation, we presented results from a theoretical approach that combines the leading 
order perturbative QCD baseline calculation of inclusive light hadron and photon production with cold nuclear matter 
effects that arise from the elastic, inelastic and coherent parton scattering in large nuclei. Our numerical simulations 
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FIG. 4. Predictions for the nuclear modification factor R p pb as a function of transverse momentum pr for photon production 
in minimum bias p+Pb collisions at the LHC energy of ^/s — 5 TeV for rapidities y = 0, y = 2, and y — 4. 
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included the isospin effect, Cronin effect, cold nuclear matter energy loss and resummed QCD power corrections to 
the leading twist results. We found that this approach can describe quite well the nuclear modification factor in 
\/s = 200 GeV d+Au collisions at the RHIC for both photon and light hadron production, including mid and forward 
rapidities. Detailed theoretical model predictions were then presented for inclusive particle production relevant to the 
recently completed p+Pb run at y/s — 5 TeV at the LHC. Our results for the nuclear modification of light hadrons and 
photons at forward rapidities fall between the ones based on EPS09 shadowing parameterization and CGC approaches. 
These predictions will soon be confronted by new experimental data to help constrain the magnitude and clarify the 
origin of cold nuclear matter effects. Within the framework of our approach, such comparison will also shed light on 
the transport properties of cold nuclear matter. 
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